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Рассмотрены результаты разработки математической модели первичных двигателей, позволяющей без декомпозиции с высо- 
кой достоверностью моделировать все виды и типы первичных двигателей с детальным учетом разнообразных систем регули- 
рования котлоагрегатами и турбинами. Приведены сведения об апробации и практическом использовании результатов. 


Одним из наиболее сложных элементов совре- 
менных электроэнергетических систем (ЭЭС) яв- 
ляется первичный двигатель (ПД) электрических 
синхронных генераторов, которым обычно служит 
гидравлическая или паровая турбина со своими си- 
стемами регулирования. В последнем случае в со- 
став ПД входит также парогенератор с присущими 
ему системами регулирования и вспомогательным 
оборудованием (котлоагрегат). 

Известная специфика энергосистем сводит к 
минимуму возможность получения информации о 
процессах в оборудовании и ЭЭС натурным путем, 
а их чрезвычайная сложность значительно ограни- 
чивает применение физического моделирования. В 
результате основным способом получения инфор- 
мации, необходимой для проектирования, эксплу- 
атации и развития ЭЭС, оказывается математиче- 
ское моделирование. 

Динамика первичных двигателей существенно 
влияет на нормальные и аварийные процессы произ- 
водства, распределения и потребления электроэнер- 
гии. Поэтому требование к уровню адекватности ма- 
тематической модели ПД постоянно возрастает. 

Применяемые в различных программах расчета 
процессов в ЭЭС математические модели первич- 
ных двигателей весьма упрощены и отражают не- 
кий обобщенный ПД. В результате достоверность 
информации, полученной при моделировании, 
оказывается низкой. Использование такой инфор- 
мации при проектировании и эксплуатации энерго- 
блоков и ЭЭС в целом прямо или косвенно стано- 
вится причиной неправильной настройки систем 
регулирования, противоаварийной автоматики и 
нередко приводит к тяжелым авариям. В связи с вы- 
шеизложенным поставлена задача разработки мате- 
матической модели ПД, позволяющей достаточно 
точно моделировать любой его конкретный тип. 

Решение поставленной задачи осложнено тем, 
что в настоящее время в ЭЭС находится в эксплуа- 
тации множество различных типов турбин и котло- 
агрегатов с разнообразными системами и законами 
регулирования. Хотя передаточные функции от- 
дельных элементов и звеньев этого оборудования 
известны, разработка цельной математической мо- 
дели для конкретного типа оборудования и компо- 
новки ПД оказывается уникальной. Это подтвер- 
ждается многочисленными исследованиями осо- 
бенностей и специфики работы парогенераторов с 


прямоточными и барабанными котлоагрегатами, 
применяемых систем и законов регулирования эт- 
ими агрегатами, конденсационных турбин, в том 
числе с промежуточным пароперегревом, турбин с 
противодавлением, турбин с промышленными и те- 
плофикационными отборами пара, гидротурбин и 
разнообразных систем регулирования перечислен- 
ными турбинами [1—5], которые показывают нали- 
чие у них как общих, так и принципиально отли- 
чающихся свойств. Таким образом, математическая 
модель даже самого сложного первичного двигате- 
ля не включает в себя более простые, и ее разработ- 
ка не решает задачи в целом. Следовательно, встает 
задача разработки математической модели для каж- 
дого существенно отличающегося типа ПД. При 
этом необходимо учитывать для каждого вида пер- 
вичного двигателя тенденции их модернизации, в 
том числе возможность значительных изменений. 

Множество сложных математических моделей, 
в которых, не смотря на существенные различия, 
есть и общие фрагменты или сравнительно просто 
адаптируемые, является не оптимальным и не 
удобным для практического использования. Поэ- 
тому возникает еще и задача синтеза, на основе 
разработанных для каждого вида и типа ПД мате- 
матических моделей, универсальной модели. 

Опуская промежуточные результаты анализа, 
разработки и синтеза ниже приводится конечный 
результат решения поставленной задачи. Наиболее 
наглядной и исходной формой сложной динамиче- 
ской модели для всех возможных методов, способов 
и средств расчета является операторная схема, кото- 
рая не только иллюстрирует математическое содер- 
жание модели, но и показывает структуру моделиру- 
емых объектов, что позволяет легко ориентировать- 
ся в ней и адаптировать ее для конкретных целей. 
Такая операторная схема разработанной математи- 
ческой модели первичных двигателей представлена 
на рис. 1, где соответствующими передаточными 
функциями отражены: К ш - загрузка энергоблока, 
определяющая в регуляторе мощности (РМ) задава- 
емую энергоблоку мощность А № К шу - уровень по- 
слеаварийной разгрузки (догрузки) энергоблока 
іАРпау, воздействующей через РМ и через электро- 
гидравлический преобразователь (ЭГП); 

Рэ _ ^ІІМ 

Рт 1 + Т имР 

- измеритель мощности синхронного генератора (СГ); 
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СО и 


К \ 


ЧК 

®ДЧ 1 + Т ЧК Р 


Ад 


Г ЗД Р 

1 +о 


- звено динамической коррекции МдК РМ; А Др - ка- 
нал коррекции МдК РМ по давлению свежего пара; 

Р т . К т , 


со,, 


со 


*ич 

1 + Т' ш р + Т" 1Ч р- 


- измеритель скорости вращения турбины; кор- 
рекция БдК РМ; 


со 


ИЧ 

СО 


к ш 

1 + Т 'тР + Т шР 2 


ИЧ-Г 1 х ич/ 

- измеритель скорости вращения турбины; 


со,. 


Т ічкР 


где: АР'дир=ѵаг — часть АР Ш? , снимаемая по экспо- 
ненциальному закону г АИР =ѵаг; 


- частотный корректор медленнодействующего 
контура (МдК) РМ; 


Р К 

1 ЭГП _ /Ѵ ЭГП 


1 + т эш р 


- общий канал БдК РМ; 


Р ^ 1 + Т?мР 

- общий канал МдК РМ; К ш - общий канал регу- 
лирования мощности конденсационных турбин в 
режимах на скользящем давлении пара и регулиро- 
вания давления пара перед турбиной «до себя», 
включая стерегущий режим, а также регулирования 
мощности турбин с противодавлением, турбин с 
промышленными и теплофикационными отборами 
пара; _Кдр ЧМ - многофункциональное управление, в 
том числе для воспроизведения системного АРЧМ; 
р 1 


К и МУ1 = 

“-МУТ 11 п , , гг 

Р&М 1 + Т мут р 

- механизм управления турбиной; К п и А нкн - ка- 
нал начальной динамической коррекции неравно- 
мерности быстродействующего контура (БдК) РМ, 
посредством которого через ЭГП формируется 
корректирующее воздействие РНК с ограничени- 
ем, согласно уравнению /нк = (/э _ ^пРпп)^жн> г Д е: 
Р э - мощность СГ, а р пп - давление пара за проме- 
жуточным пароперегревателем; А скн - статическая 
коррекция неравномерности БдК РМ; 

_ 1 + Г дм Р 

р, 1+0 

- звено динамической коррекции БдК РМ; К п А - 
динамическая коррекция БдК РМ; 


^РС=- 


- регулятор скорости (РС), где: о=ѵат - статизм; 
К„, - задание уставки щ РС; 

Рп 1 
Р, 1 + Т Ц Р 

- промежуточный золотник регулятора турбины; 
+К т и ±А РДН - каналы акцентного управления регу- 
лирующими клапанами (РК) части высокого давле- 
ния (ЧВД) и части низкого давления (ЧНД) турбины 
с промышленными и теплофикационными отборами 
пара при осуществлении связанного регулирования, 
обеспечивающего автономность управления скоро- 
стью вращения турбины и давлением; К щ и К щ - ко- 
эффициенты передачи каналов акцентного управле- 
ния РК ЧВД и ЧНД турбины по скорости вращения; 

Ч.,.о _ 1 

Чзво 1 + 7 ™ Р 

- открытие окон золотника сервомотора РК ЧВД; 

Чво _ 1 

Рщ, Т »оР 

- перемещение поршня сервомотора РК ЧВД на 
открытие; 


Р я 


1 


Чзвз 1 + Т жР 

- закрытие окон золотника сервомотора РК ЧВД; 


Рп 


1 


Р, вз Р т Р 

- перемещение поршня сервомотора РК ЧВД на 
закрытие; 

Р, со _ 1 

Чзсо 1 + 7, зсо Р 

- открытие окон золотника сервомотора РК части 
среднего давления (ЧСД) или ЧНД, или их обоб- 
щения, которое обычно также обозначают ЧНД; 


®ЛЧ 1 + РщкР 

- динамическая частотная коррекция БдК РМ; 
К ? - многофункциональное управление разгрузкой 
через БдК РМ, в том числе противоаварийной ав- 
томатической импульсной разгрузкой (АИР), осу- 
ществляемой путем функционального изменения 
Ар, обеспечивающего формирование: 

I 

Д Р = др(°) . р Ты№ 

ш АИР ш АИР с > 


Ра 

Р, с 


1 


т Р 

СО г' 


- перемещение поршня сервомотора РК ЧСД, 
ЧНД на открытие; 


/А 

Чз< 


1 


і + ^зсзЧ 


- закрытие окон золотника сервомотора РК ЧСД, 
ЧНД; 
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АА, 
/К . 


1 


Т Р 

СЗ-Г 


- перемещение поршня сервомотора РК ЧСД, 
ЧНД на закрытие; К ожв и ^ ожс - коэффициенты же- 
сткой отрицательной обратной связи систем упра- 
вления РК ЧВД и ЧСД, ЧНД, соответственно; 

_ Т шР 
И, 1 + Т'„р 

- гибкая отрицательная обратная связь системы 
управления задвижкой направляющего аппарата ги- 
дротурбины, функции которой в этом случае возла- 
гаются на систему управления РК ЧВД с соответ- 
ствующими изменениями параметров передаточных 
функций, а аналогичная система ЧСД, ЧНД при 
этом не используется; - степень открытия РКЧВД 
паровой турбины или задвижки направляющего ап- 
парата гидротурбины; +К ЮВ - многофункциональное 
задание предварительного положения РК ЧВД, за- 
движки направляющего аппарата или различного 
рода возмущающих воздействий типа дребезга и др.; 
С с - степень открытия РКЧСД, ЧНД паровой турби- 
ны; ±Д,, д - многофункциональное задание предва- 
рительного положения РКЧСД, ЧНД или различно- 
го рода возмущающих воздействий и др.; 

РпП _ 1 

Дпп Т пп р 


- промежуточный пароперегреватель, где: Д пп - 
расход пара через пароперегреватель; Д в - расход 
пара через ЧВД, равный произведению \ в на соот- 
ветствующее давление пара, в зависимости от типа 
турбины и режима работы согласно функциональ- 
ной операторной схеме, или расход воды при моде- 
лировании гидротурбины, равный произведению 
<^ в (степень откры тия задвижки направляющего ап- 
парата) на л/Т +Н В в соответствии с законом истече- 
ния жидкости через задвижку, где: Н в - относитель- 
ное превышение напора воды перед задвижкой на- 
правляющего аппарата, определяемое с учетом яв- 
ления гидроудара (ГУ) передаточной функцией [3] 


Я 


Д 


= -24* 


1 - е~ ТтР 
1 + е~ ТгР ’ 


ш аѴ 

где. 4* = — — - ударная характеристика водовода, 
2Н 8 

а - скорость распространения волны ГУ, V - ско- 
рость потока воды, Н е - номинальный напор воды, 

2 Ь 

т г = 

а 

- фаза ГУ для трубопровода длиной Ь\ 

Д вд _ 1 

Дв І+Твд Р 

- паровой объем ЧВД, где: Д вд - расход пара через 
этот объём; К щ - доля мощности и соответственно 
момента турбины (М т ) за счет ЧВД; К а - часть до- 
ли мощности и соответственно М т за счет ЧСД; 


Дсд_ 1 
Дс 1 + 7> 

- эквивалентный паровой объема ЧСД, располо- 
женный за контуром промперегрева, где: Д сд - рас- 
ход пара через этот паровой объем, Д с - расход па- 
ра через ЧСД, равный произведению ос на соответ- 
ствующее давление пара согласно функциональ- 
ной схеме; К в - доля мощности и соответственно 
М т за счет ступеней ЧСД, расположенных за ука- 
занным выше эквивалентным паровым объемом; 

Дне _ 1 

Дсд 1 + 7 Н сР 

- эквивалентный паровой объем, расположенный 
за РК ЧНД, где: Д нс - расход пара через ЧНД; К ш - 
доля мощности и соответственно М т за счет ЧНД; 
К Мт - многоцелевое управление, в частности для мо- 
делирования различного рода возмущений, а также 
обнуления момента гидротурбины М г =К Мг Д в (1+Д,), 
при моделировании паровой турбины; К ш - вос- 
произведение разнообразных аддитивных возмуще- 
ний М т ; К п , К ру , Кр, _К, Н , К Чв , К к од, Кцр, Кд, Кд в , К ш — 
многофункциональные коэффициенты, позволяю- 
щие, в частности, задавать тип и режим работы мо- 
делируемого первичного двигателя: 

Рк _ 1 
бдк Т кР 

- парогенератор (котлоагрегат (КА)), где: 0 ДК - ко- 
личество тепла, необходимое для производства 
расходуемого количества пара Д д с давлением р к ; 

О к 

УДъ = е ~ т гР 1 — _ тракт топливоподачи, где: 

Дд І + Т'тР 


Дд 


е 


-*тР. 


бдт 


К т 

топка; 

1 + Т л р 


Рдк _ ^дкР 
Р 1 + ПкР 


- динамическая коррекция контура регулирования 
производимого КА тепла по давлению; К т? - про- 
порциональная часть пропорционально-инте- 
грального (ПИ) регулятора топлива (главного регу- 
лятора) КА; 

Дди ^ІІР 

Ддр т ѵ?Р 

- интегральная часть ПИ - регулятора топлива КА; 

Дсн _ ^сн 
Од 1 + ГсиР 

- коррекция производительности КА из-за изме- 
нения режима работы оборудования собственных 
нужд, обусловленного отклонениями частоты ю д ; 
Кд, - задание паропроизводительности Д, д КА и ее 
изменения Д Дз ; АГ РМД - канал РМ управления произ- 
водительностью КА; К в і — перепад давления пара 
перед турбиной, определяемый уравнением 
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Рис. 2. Процесс противоаварийной разгрузки энергоблока 800 МВт Сургутской ГРЭС- 2, где: М т ~ момент турбины, Р г - актив- 
ная мощность генератора, ЭеНа ~ взаимный угол генератора 

аварийной разгрузкой энергоблоков, что позволяет 
более достоверно оценивать динамическую устой- 
чивость ЭЭС и оптимально настраивать средства 
противоаварийной автоматики. 

Адекватность рассмотренного моделирования 
различных видов и типов первичных двигателей под- 
тверждены результатами его использования в составе 
созданных в НИЛ «Моделирование ЭЭС» гибрид- 
ных моделирующих комплексов, в частности для 
ОАО «Тюменьэнерго», на котором проведены ус- 
пешные исследования противоаварийной разгрузки 
энергоблоков Сургутских ГРЭС-1 и ГРЭС-2 [6]. 
Фрагмент этих исследований иллюстрирует рис. 2. 

Следует также отметить большое значение рас- 
смотренного моделирования первичных двигате- 
лей для создания надежных и эффективных отла- 
дочно-тренажерных комплексов АСУТП энерго- 
блоков [7]. 

Выводы 

1. Разработанная математическая модель первич- 
ных двигателей синхронных генераторов позво- 
ляет адаптировать ее для моделирования всех ви- 
дов и типов используемых первичных двигателей 
с учетом их систем и законов регулирования. 

2. В полученной математической моделе первич- 
ных двигателей отображены все значимые эл- 
ементы оборудования, что служит объективной 
гарантией повышения уровня адекватности мо- 
делирования. 

3. Основные результаты разработки подтвержде- 
ны практикой их использования при эксплуата- 
ции Гибридного моделирующего комплекса 
Тюменской энергосистемы. 


Дрт = А^Д 2 по сравнению с давлением на выходе КА 
р к , где: Д - расход пара через турбину; 

Рт _ К кт 

Ркт \ + Т ш р 

- паропровод между КА и турбиной; К я - обоб- 
щенный измеритель давления пара перед турби- 
ной; К р з - задание давления пара перед турбиной; 
К л , - обобщенный измеритель расхода пара через 
турбину; 

Ддц _ ^АД 
Ддп 1 + 7 ед Р 

- изменение потребления пара из коллектора 
вследствие каких-либо возмущений у потребителя 
Дагъ К птр - взаимосвязь между изменением давле- 
ния в коллекторе А р РВК и отклонением расхода па- 
ра АД; К рг - задание уставки давления в коллекторе 
Рк»; Кш — общий канал главного ПИ-регулятора 
КА; К РВ - пропорциональная часть ПИ-регулятора 
давления КА; 

Ррр _ Кррк 

Р&с Т рвк р 

- интегральная часть ПИ-регулятора давления КА. 

Синтезированная модель учитывает практиче- 
ски все основные процессы в первичных двига- 
телях, в результате чего значительно повышается 
точность моделирования динамического баланса 
генерируемой и потребляемой мощности в ЭЭС, а, 
следовательно, процессов изменения частоты и пе- 
рераспределения генерации между энергоблоками 
и электростанциями с учетом их колебательной 
устойчивости. Значительно точнее воспроизводят- 
ся динамические переходы, связанные с противо- 
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